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Abstract 

We present analytical and numerical results for the specific heat and 
susceptibility amplitude ratios in parallel plate geometries. The results 
are derived using field-theoretic techniques suitable to describe the sys- 
tem in the bulk limit, i.e., (L/£±) 3> f, where L is the distance between 
the plates and £± is the correlation length above (+) and below (-) the 
bulk critical temperature. Advantages and drawbacks of our method are 
discussed on the light of other approaches previously reported in the lit- 
erature. 

PACS numbers: 64,10+h, 64.60.Fr, 64.60.Ak 
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I. INTRODUCTION 



Since the advent of modern scaling concepts and renormalization- 
group techniques the study of finite-size and surface effects on the be- 
havior of systems near or at criticality has attracted the attention of a 
number of investigators 1 . 

Fixing our interest in the case of a system confined between two in- 
finite {d — l)-dimensional parallel plates distant L from each other, we 
may classify three well defined distinct regions in this problem. The first 
one, where the scaling variable (L/£±) ^> 1, is characterized by the dom- 
inance of bulk over surface and finite-size effects and the physics is quasi 
d-dimensional. Here, £± specifies the critical correlation length above 
(+) and below (— ) the bulk critical temperature T c . The second region, 
where (L/£±) <C 1, the system behaves as a quasi (d — l)-dimensional 
object. Finally, for (L/£±) ~ 1, the physics interpolates between a quasi 
d- and a quasi (d — l)-dimensional system. A full description of the 
system should therefore unveil very interesting crossover behaviors. 

According to the region and phenomenon of interest, different field- 
theoretic techniques have been devised to deal with such systems. In- 
deed, Diehl and Dietrich 2,3 successfully implemented these techniques to 
study critical and multicritical phenomena near surfaces within a finite 
momentum cutoff regularization scheme. The use of dimensional reg- 
ularization was shown 4 ' 5 to simplify the computational procedure and 
allowed the study of ordinary 2 and special transitions 3 through standard 
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4 -field theories under Dirichlet (DBC) and Neumann (NBC) boundary 
conditions, respectively. The former mimics very strong repulsive forces 
at the surface, thus preventing order at it (a parameter c, which mea- 
sures these forces 2 , has fixed point value c* = oo), whereas under the 
later boundary condition both the surface and the bulk go critical simul- 
taneously. The special transition is in fact a multicritical point 3 , c* = 0, 
where the two lines describing systems with repulsive (c > 0) and at- 
tractive forces (c < 0) at the surface, meet. In the later case, namely the 
extraordinary transition 5 ' 6 , the surface undergoes a second order tran- 
sition before criticality sets in the bulk. Moreover, it has been shown 7 
that a scaling description holds so that the critical exponents associated 
with excess surface singularities may be expressed completely in terms 
of bulk exponents. However, it has also been shown 7 that fluctuations 
may induce divergences at the surface and in these cases local quantities 
and associated exponents must be defined resulting in new scaling rela- 
tions. Therefore, in treating these quantities at the multicritical point, 
Neumann boundary conditions are valid only at the mean-field level 7 . 

On the other hand, in order to properly describe finite size effects 
using field-theoretic techniques in critical systems subject, for example, 
to periodic boundary conditions (PBC), Brezin and Zinn- Justin 8 and, 
independently, Rudnick et al 9 , introduced a method in which the zero- 
momentum component is isolated whereas the other non-zero modes are 
treated perturbatively. This methods has been largely used 10 ' 11 ' 12 and 
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generalized to study different boundary conditions. More recently , 
some difficulties to treat critical systems below T c using this technique 
have been circumvented. 

Both finite size and surface effects are simultaneously present, except 
in special circumstances such as for PBC where surface effects do not 
contribute. At T = T c , where Casimir forces manifest, these contribu- 
tions compete in a very special way, and powerful tools and methods 
such as conformal invariance 14 and elaborated perturbation techniques 15 
have been used to study this regime and the approach to J^ 12 ' 13 ' 15 ' 16 . 

In this work we shall calculate specific heat and susceptibility ampli- 
tude ratios using field-theoretic and G-expansion methods 16 particularly 
suitable to describe systems in the first regime mentioned above in which 
bulk behavior dominates over surface and finite size contributions. The 
reported results complement previous studies 16 and shed some light on 
the approach to bulk criticality as the distance L between the plates 
increases. 

For a system of volume V = AL, where A is a [d — l)-dimensional 
surface (layered geometry), the following asymptotic scaling form for the 
singular part of the free energy density holds 17 : 

/(N,L)~^F(L|t|7a)=^|tr + y s ^^ + (5/(L|t|7a) , (1) 

where t = (T — T c )/T c , v is the bulk correlation length exponent and 
a is the only non-universal metric factor. Using the hyperscaling rela- 
tion dv = 2 — a, one identifies the first term (proportional to y^) as 
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the bulk contribution, the second one (proportional to y s ) as the excess 
surface term and the last one as a finite-size correction term. In the 
limit (£±/L) <C 1 one expects exponentially small corrections from 8f, 
whereas for (£±/L) ^> 1 it compensates the bulk and surface contribu- 
tions and gives rise to the Casimir effect 13-15 at T = T c . 

From the above scaling assumptions the specific heat and susceptibil- 
ity should behave as 

C(t,L) ~ \t\- a A±(L\t\ u /a) , (2) 
X (t,L)~\trC ± (L\t\»/a) , (3) 

where a and 7 are the bulk critical exponents, but even in the regime 
L\t\ p >> 1 one expects that excess surface and finite-size contributions 
modify their critical amplitudes in a non-trivial manner. In fact, the 
ratio of these amplitudes are quite sensitive in identifying the universality 
class of a critical system, particularly in numerical simulations 18 where 
one has to control both corrections to scaling and surface and finite size 
effects. 

In Section II we explain our method and derive both the renormalized 
free energy and the equation of state from which the above quantities 
can be calculated. Finally, in Section III a discussion of the results and 
conclusions are presented. 
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II. SPECIFIC HEAT AND SUSCEPTIBILITY CRITICAL AMPLITUDES 



In this section we shall use field-theoretic and renormalization-group 
techniques to calculate the amplitude ratios of C and x m layered geome- 
tries. We shall keep close contact with standard bulk 4 — field theory 19 
and whenever necessary to deal with the finite size of the system we em- 
ploy methods 16 which are particularly suitable in the regime (L/£±) ^> 1. 

A. Renormalized Free Energy and Boundary Conditions 

We start by writing the expression for the one-loop renormalized 
Helmholtz free energy density at the fixed point associated with the 
bulk critical behavior of the system: 

F(t, M; L) = -tM 2 + -u*M A + - (t 2 + uHM 2 + -u* 2 M A \ I sp 
2 4! 4 V 4 J y 

+^ E / d d - l q In [1 + (l/2)u*M 2 /(q 2 + k 2 + t)\ . (4) 

In the equation above t, M (to = Z^t, Mo = Z l J 2 M^j are the renormal- 
ized (bare) reduced temperature and order parameter, respectively, Z^, 
Z^ are renormalization functions, u* is the dimensionless renormalized 
coupling constant of a continuous (bulk) 4 theory at the fixed point, <f is 
a (d— l)-dimensional wave vector along the direction parallel to the plate, 
Kj = Tij/L are the eigenvalues of the kinetic energy operator satisfying 
proper boundary conditions (see below) and I sp =G _1 [1 + (G /2)]+0(g) 
is the one-loop integral of a bulk 4 theory evaluated at the symmetry 
point using dimensional regularization, where G= 4 — d. Notice that 



taking the limL — > oo in Eq.(4) one obtains the standard expression for 
the d-dimensional 4 one-loop renormalized free energy. 

In deriving F(t,M; L) we have considered that the local field of 
a (j) A theory may satisfy periodic (PBC), Neumann (NBC) or Dirich- 
let (DBC) boundary conditions, defined by: (/>(p,z) = 4>{p, z + L), 
(d/dz\ z= o)<f>(p,z) = (d/dz\ z=L )(f>(p,z) = and <f>(p,z = 0) = <f>(p,z = 
L) = 0, respectively, where p is a [d — l)-dimensional position vector 
perpendicular to the z(dth) direction. It then follows that the sum in 
Eq.(4) have values j = 0, ±1, ±2,..., for PBC, j = 0, 1, 2, for NBC, 
and j = 1,2, for DBC, respectively. The local field is Fourier trans- 
formed in the form 16 

</>(p, *) = E(27r) 1 - d / d d - l q exp{iq.p)^{q) U] {z) , (5) 
j 

where 4>j(q) are plane waves parallel to the plate and Uj(z) are eigenfunc- 
tions of the kinetic energy operator (—d 2 /dz 2 ) with eigenvalues k|. The 
bare order parameter Mq is thus the expectation value of the local field 
above. We call attention that the usual counterterms of a bulk 4 theory 
are used to renormalize the free energy and that the boundary conditions 
are implemented on the bare vertex functions. Details of the Feynman 
rules involving propagators and vertices can be found in Ref.(16). 

B. Specific Heat Amplitude Ratio 

Since the vertex function T^ ' 2 ) is additively renormalized, the criti- 
cal behavior (singular part) of the specific heat is calculated using the 
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expression : 



V - -(0,2) 



C = A\t\- a = — B{u*) - T% > , (6) 



a 

where B(u*) is the inhomogeneous term of the renormalization group 
equation for T^' 2 ^ and 

rf 2) = |^(«>M;L) . ( 7 ) 
For T > T c , M = 0, and we find, using Eqs.(4) and (7) 

rg' 2) (T>r c ) = -±E/-i^L- 5 + h , (8) 

2L J J (q 2 + *) + £) 2 

whereas for T < T c we use the value of M at the coexistence curve, 
namely u*M 2 = —6t, to obtain 

rr>(r<r c) = -l-fW 2 +2i sp . (9) 

L j J (q 2 + K 2 + 2\t\) 



u* 



The one-loop integrals are evaluated using dimensional regularization 
and some useful formulae 16 ' 20 to sum infinite series. We thus obtain for 
the boundary conditions of interest: 



Lt^' 2 ' 

2 a TT 



Lj J (q 2 + K 2 + t) 2 eV 2) 

+27T- 1 / 2 (-) r(d/2)r[(5 - d)/2)]s*n[7r(5 - d)/2]f {5 _ d)/2 

IT 1 / 2 

+r_r(d/2)r[(5 - d)/2] p (d - 5)/2 /^] , (10) 

where 5^/ (27r) d = 1, 5^ being the area of the d-dimensional unity sphere, 
t = t+ (l/2)w*M 2 = tit = 2\t\) for T > T C {T < T c ), a = 1 for PBC, 



(7 = for both NBC and DBC, r = 0, +1, -1 for PBC, NBC and DBC, 
respectively, and 

, x too (u 2 — a 2 ) a du , s 

fa(a) = / V ft) \ - , a < 1 • 11 

exp(Z7TU) — 1 

Now using the G-expansion 19 for the non-singular part of the specific 

heat, -(v/a)B{u*) = (3/2 g) + (295/108) +0(g), we find the amplitudes 

above and below T c : 
3 



^4+ = 



A_ = 



G 2 
6 

e~2^ 



47 2 2 - a 

7 2 2 - 
! - G 54 +G 3 



, (Lt l l 2 \ T 7T 

/l/2 v Wj + G 3L^ 



2 2 -^ /V2^|t| 1/2 \ 
~3~ /l/2 1, 2-tt J 



+ 0(G 2 ) , (12) 
+ 0(e 2 ) , 



T7T 



Z^/2L\t\ l l 2 



(13) 



where a =G /6 + C(g 2 ). 

From the above equations we finally obtain the specific heat amplitude 
ratio: 

A oa 



A 2 a 

-^ = -[i+ef( x )+eS A (x)\ + 0(e- 



(14) 



where f(x) and Sa (x) are given by 

r>2-<7 

/(*) = ^r- 



/i /2 (x/2%) - / 1/2 (V2i/2"7r)J , (15) 

^^^t 1 "^) ' (16) 
and x = L/£, with £ = |^| -1 ^ 2 . For x — > oo, we can use the asymptotic 

limit 16 ' 20 for f 1/2(0) in Eq.(ll) and write f(x) in the more simplified form 



>1-(T 



(27T) 



1/2 



exp(— 2 



l-cr 



exp(- v / 2 2 1 ~ a 



x 



(2 1 --x) 1 / 2 ( v / 2 2 1 --x) 1 /2 



00 



(17) 
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C. Susceptibility Amplitude Ratio 



Using Eq.(l) we obtain the following renormalized equation of state: 



dF 1 1 / 1 N 

H R = — = tM + -u*M 3 + -u*M t + -u*M 2 
K dM 6 2 V 2 



x 



1 



d d - l q 



Lp ^ U (f + «*) (g 2 + + t + ± W *M 2 ) 



(18) 



The one-loop integral is calculated similarly as for the specific heat: 



d d - 



q 



, 2 ^ 2 w 2 ^ 2 — TV = ^i 2 - (e /2)]r(e /2) jf 1 ^(te)- e / 2 

+27T-V2 4 T(d/2)r[(5 - d)/2]sm[7r(5 - d)/2] 



L 



X I dx f($-d)/2 



L{tx) 



(5-d)/2' 



2 ct tt 



+r^r(d/2)r[(5 - cO/2] / rfx [(te)^ 5 )/ 2 /L 



(19) 



By noticing that the first term in the right-hand side of Eq.(19) may 
be written as [g -1 — {1/2) int ], and using the G-expansion representation 
for I sp , we obtain, to first order in G, 

H R = tM + ^u*M 3 + ^u*Mt{{! + int) 



L{yt) 



1/2 



2 ff 7T 



Tr r 1 / 2 
L 



(20) 



where w* = (2/3) G +C(g 2 ) and f = £ + (l/2)w*M 2 . 

The susceptibility amplitudes are then readily calculated from 



v -i _ (ritnV 1 - r m - ^ 
x -{c\t\ ) -i R - m 



(21) 
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As before, for T > T c , M = 0, and for T < T c we use u*M 2 = -6t. The 
amplitudes above and below T c thus read: 

,1/2^ 



C+ = 1 -| + |/ 1 ^ /l /2 

C_ = i{l-|(4 + £n2) + f/ o 1 ^/ 1 



T7T 



+ G— + (9(G 2 ) 
7T / OX 



(22) 



x(2y) 



1/2' 



2 CT 7T 



+ G 



T7T 



6^ 



+ C(G 2 ) 



(23) 



where x = L/£. 

Using Eq.(ll) and performing the integrations in y, we find the sus- 
ceptibility amplitude ratio: 



^ = e h(x)+ e S c (x) + G(e 2 ) , 

where 7 = 1 + (<E /6) + C(g 2 ), /? = 1/2 - (g /6) + C(g 2 ) and 



(24) 



h{x) = 



2" +1 7r 
3x 



24 



— / dug{a/u) + v2 / dug(b/a) 



with a = V2x/2 a 7i, b = x/2 <7 n and 

wcosfarc sm(c/w)] 



(25) 
(26) 



fl(c/«) = 



(27) 



III. DISCUSSION AND CONCLUSIONS 

First, we should point out that the main step in our approach is the 
representation 16,20 used to evaluate the discrete sums in Eqs.(10) and 
(19). It has proved very useful in different field-theoretic contexts 20 and 
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here it clearly helps to split the bulk, surface and finite size contributions, 
as required by scaling [see Eqs.(l-3)], in a rather simple way, though its 
range of effectivenes precludes direct access to the Casimir effect. 

Second, our starting renormalized free energy, Eq.(4), does not con- 
sider any distortion of the order parameter profile, i.e., our descrip- 
tion is restricted to calculating the effect of the boundary conditions 
on bulk quantities as a result of fluctuations, i.e., in the amplitude 
ratios, Eqs. (14) and (24), excess surface and finite size contribu- 
tions are of 0(g). Nevertheless, we observe that if the excess sur- 
face contributions for NBC in Eqs. (12) and (13) are isolated, we find 
(A + /A_) s = 2~ 3 / 2 + 0(e), which is the same result derived in Ref.(5) 
for the special transition (the excess surface specific heat exponent is 
a,s = a + v). This is so because in this particular case there is no 
distortion of the order parameter at the mean-field level. Notice also 
that, above T c , (A + ) sp /(A + ) ord = -1 + 0(e 2 ), where {A + ) sPj0rd refer to 
the specific heat amplitudes at the special (NBC) and ordinary (DBC) 
transitions, a result already derived in Ref.(4). As for the excess surface 
contributions for the susceptibility amplitudes we also notice that the 
last term in Eq.(20) is consistent with the result found in Ref.(7) for the 
special transition above T c (75- = 7 + 1/), since again no distortion of the 
order parameter is necessary in this case. 

From the discussion above and the derived results in Section II, par- 
ticularly Eqs. (14-17) and (24-27), it is clear that in the regime L/£± ^> 1, 
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and to first-order in an 6-expansion, the specific heat and suscepti- 
bility display singularities well described by bulk exponents, but with 
amplitudes sensitive to the boundary conditions which manifest as ex- 
cess surface and finite size contributions. Notice also that these fluc- 
tuation effects result quite effectively from the difference between the 
amplitudes of the correlation length above and below T c , which satisfy 

(&,+/&,-) = >/2 + 0(e). 

In Figs. 1, 2 and 3 we plot the scaling functions f(x) and h(x), 
numerically evaluated in d = 3, as defined by Eq.(ll, 15, 25-27). They 
both decay very rapidly to zero as x = Lt 1 ^ 2 increases, in agreement with 
the asymptotic result 16 ' 20 for fi/2(x), x — > oo. In fact, for PBC and x = 7 
we find f p (x) ~ 2, 7 x 1CT 4 both by using the numerical estimate and the 
asymptotic result predicted in Eq.(17). For this value of x it is indeed 
expected 18 that these corrections to the bulk limit are indeed negligible. 
Notice that for NBC the magnitude of the scaling functions are the same 
in our one-loop approximation, in agreement with Ref.(16). However, a 
two-loop calculation shows 16 that they slightly differ if the same regime 
of validity applies. However, as x — > our approach does not correctly 
describe the Casimir effect, as shown in Figs. 1 and 2: f(x) diverges and 
h(x) approaches zero, whereas in a correct treatment 13-15 both tend to 
a constant value, the Casimir amplitudes for each case. This failure for 
x < 1 has already been pointed out by Nemirovsky and Freed 16 . Here, 
our results deal with issue in a more quantitative way, thus evidencing 
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the advantages, limitations and drawbacks of the method. 

For comparison, we also plot in Fig. 1 the excess surface contribution 
for the cases of both NBC and DBC, i.e., \S A (x)\ = \S c (x)\ = S(x). 
As clearly seen from Figs. 1 and 2, as x — > oo these contributions are 
the leading ones 11 ' 15 modifying the amplitudes of both the bulk specific 
heat and susceptibility. Notice that they having the same magnitude 
and decay as x~ x ( this is fortuitously true only because these effects are 
treated as fluctuation contributions to the bulk limit, see Eqs. (16,26) ). 

Finally, in Fig. 4 we plot the difference between the excess surface 
contribution and the scaling function for NBC. We see that this difference 
"almost" saturates as x — > 0, as expected in the Casimir effect, but lostly 
it diverges for very small values of x (a x^ 1 dependence is in fact expected 
from Eq. (11) as x — > 0, but an extra \nx contribution precludes a good 
description of the Casimir effect). 

In summary, we have presented a field-theoretic description of the 
approach to bulk criticality in parallel plate geometries, in which excess 
surface and finite-size contributions appear as a result of fluctuations 
and depend on the boundary condition imposed on the system. Despite 
the fact that other more general methods to deal with finite systems do 
exist, our approach is probably the simplest one in the regime L/£± ^> 1. 
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FIGURE CAPTIONS 



Fig. 1. Amplitude of the excess surface contribution (solid line), from 
Eqs. (16) and (26), and scaling functions f p (dashed line) and f^,D 
(dotted line) for periodic, Neumann and Dirichlet boundary conditions, 
respectively, numerically evaluated using Eqs. (15) and (11) in d = 3, 
as functions of x = L/£. 

Fig. 2. Same as in Fig. 1 for 2 < x < 7. 

Fig. 3. Scaling function h p (dashed line) and Hn,d (dotted line) for 
periodic, Neumann and Dirichlet boundary conditions, respectively, nu- 
merically evaluated using Eqs. (25) and (27) in d = 3, as function of 
x = L/C 

Fig. 4. Difference (apart from a minus sign) between the excess sur- 
face contribution and the scaling function fjy for Neumann boundary 
condition as function of x = L/£. 
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